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Abstract 

We consider the reversible random sequential adsorption of line segments on a one-dimensional 
lattice. Line segments of length I > 2 adsorb on the lattice with a adsorption rate K a , and leave 
with a desorption rate We calculate the coverage fraction, and steady-state jamming limits 
by a Monte Carlo method. We observe that coverage fraction and jamming limits do not follow 
mean- field results at the large K = K a /Kd » 1. Jamming limits decrease when the length of the 
line segment I increases. However, jamming limits increase monotonically when the parameter K 
increases. The distribution of two consecutive empty sites is not equivalent to the square of the 
distribution of isolated empty sites. 
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INTRODUCTION 



The irreversible adsorption of large molecules reported on systems of colloids, proteins, 



latex spheres, polymer, etc. 



QflflQQflflflflyQ 



12 



131 ] - The most simple model 



of the irreversible adsorption is the random sequential adsorption(RSA). Large molecules 
impact sequentially on the surface. If the impacting surface is empty, the molecules adsorb 
on the surface and do not detach from the surface. If the impacting surface is occupied 
by molecules, impacting molecules can not adsorb on the surface. Therefore, we expect a 
formation of a monolayer. In the RSA model the coverage fraction of the surface approaches 
a limiting value, so called, jamming limit at the long time. In the lattice model of RSA the 
coverage fraction 9(t) follows a exponential behavior as 9(t) = 9(oo) — Aexp(-Bt) where 
9(oo) is a jamming limit, A and B are constants depending on the dimensionality of the 
surface and the shape of molecules. In the continuous model of RSA or the parking-lot 
problem, the coverage fraction follows a power-law behavior as 9(t) = 9(oo) — At~ a where A 
is a constant and the exponent a depends on the dimensionality and the shape of the object. 
The RSA is a oversimplified model for the adsorption of large molecules. Indeed, there are 
a lot of effects such as the transport of molecules, the diffusion of adsorbed molecules, and 
the desorption from the surface to the bulk solution 14, 

There are many studies on effects of the desorption on the RSA for physical, chemical, 
and biological systems jl|, Q, [J. A simple example of RSA with desorption is a parking lot 
problem. Identical cars adsorb(or park) on a line(curb) at the rate K a . A certain number of 
parked cars leave a empty space that is too small to fit another car. In the irreversible model 
of the parking lot (i.e. cars park permanently on the parking lot), the density of cars reaches 

n 

to the jamming limit 9(oo) = 0.747 ■ • -[15j. In the reversible model, the cars park on a line 
at a rate K a and leave the line at a rate In the reversible model the coverage fraction 



depends on the adsorption and desorption rate. For 



arg e values of K = K a j K d the coverage 
18j |. Jin et.al. studied the adsorption- 



fraction shows two different time scales Q, 3, 
desorption process of rods on a line. They observed the logarithmic dependence of the 
coverage fraction^]. However, their approaches are based on the mean-field idea. Kolan 



et.al. reported the glassy behavior of the parking-lot model. They observed two different 
time scales of the coverage fraction by the Monte Carlo method 16| . Ghaskadvi and Dennin 
reported the reversible RSA of dimers on a triangular lattice|l2|. Krapivsky and Ben-Nairn 
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FIG. 1: The coverage fraction 0(t) versus log(i) for I = 2 with K = 10, 50, 100, 500, 100, 5000 from 
bottom to top, and K = co(solid line). 

reported mean-field results of the coverage fraction for adsorption-desorption processes [tbI]. 
They obtained the jamming limit 9{oo) ~ 1 — l/log(K), (K » 1) for the continuum 
model. The coverage fraction follows a logarithmic dependence such as 6(t) ~ 1 — l/log(t) 
for the desorption controlled limit (K a = oo,Kd = finite). In the lattice model the steady- 
state coverage fraction was given as 9 eq = l- (l/K) l / l /l (K»l) where I is the length of 
adsorbed objects. 

In this work we consider the reversible RSA process on the one- dimensional lattice by 
using Monte Carlo simulation. We observed that the coverage fraction shows three time 
scales for large values of I. However, the jamming limit at the steady-state do not follows 
the mean-field behavior for the large K. In section II we present the Monte Carlo method 
of the reversible RSA. We give results and discussions in section III and concluding remarks 
in section IV. 

REVERSIBLE RSA MODEL AND MONTE CARLO METHOD 

Consider a one-dimensional clean lattice initially. In the reversible RSA model, a line 
segment of a length I drops on the lattice at the rate K a . Adsorbed Z-mers are desorbed 
from the lattice at the rate K^. We select randomly a lattice site. If the selected site is 
empty, we try to adsorb a line segment of length I with a probability p a = K a /(K a + Kd). 
We generate a random number p(0 < p < 1). If p < p a , we try to drop a line segment on 
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FIG. 2: The coverage fraction 6(t) versus log(i) for K = 1000 with I = 2,4,8, 16,32 from top to 
bottom. 
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FIG. 3: The steady state jamming limit against the parameter log(iT) for I = 2(o) and / = 4(A). 
Solid lines represent mean field results. 

the lattice. In the trial of adsorption we check (7 — 1) consecutive neighbors. If /-consecutive 
lattice sites are all empty, Z-mers adsorb on the lattice. Otherwise, the trial is rejected and 
increase a Monte Carlo step. If the selected site is occupied by the adsorbed line segment, 
we generate a random number p. If p > p a , the occupied line segment desorbs from the 
lattice. Otherwise, the trial of desorption is rejected and increase a Monte Carlo step. 

It is not allowed to adsorb a line segment on occupied sites. So we expect a formation of 
the monolayer. Kinetics of the coverage fraction is controlled by a parameter K = K a j Kj. 
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We set K d = 1 and control the adsorption rate K a . So that, a line segment will be desorbed 
from the lattice every (K a + fQ)-trials. 

We consider a lattice of size L = 10 4 with a periodic boundary condition. The lattice is a 
one-dimensional ring with a length L. We also check the finite size effect for the large lattice 
with L = 10 5 . We calculate the coverage fraction as a function of the time. The coverage 
fraction was averaged over 100 different configurations. One Monte Carlo time corresponds 
to L adsorption/desorption trials regardless of successful or unsuccessful trials. Maximum 
Monte Carlo times were up to t = 10 5 . 

RESULTS AND DISCUSSIONS 

Consider the reversible RSA process of the line segment with the length 1 = 2. In Fig. 
1 we present the coverage fraction versus the time. For the irreversible RSA (K = oo), the 
coverage fraction saturates exponentially to a jamming limit 9(oo) = 0.868 • • ■. For a small 
K (for example K = 10) the jamming limit at the steady-state is smaller than that of the 
irreversible RSA. When the parameter K increases we observed two different linear regions 
of the coverage fraction when we plot the coverage fraction 9(t) against log(t). The coverage 
fraction increases rapidly at the early times t < t 1} and converges to the jamming coverage 
fraction of the irreversible RSA at t± < t < t 2 as shown in Fig. 2. In the former region the 
adsorption controls the kinetics of the RSA. The time t\ and t 2 increase when the length of 
the line segment / increases as shown in Fig. 2. At t 2 < t < t 3 the coverage fraction increases 
linearly up to the saturation time t 3 when we plot the coverage fraction 9{t) against log(t). 
In this region the coverage fraction is proportional to 9{t) ~ log(t). At t > £3 the coverage 
fraction saturates to the steady-state jamming limit. Steady-state jamming values increase 
when the parameter K increases. 

In Fig. 2 we gave the coverage fraction 9(t) versus log(i) for a fixed K = 1,000 and 
the different length of the line segment /. When the length of the line segment increases, 
we observe clear four stage behaviors of the coverage fraction. For small values of the 
line segment I < 8, the time t\ is very small and the coverage fraction approaches to the 
irreversible jamming value rapidly. The times t 2 and t 3 increase as the length of the line 
segment increase. Steady-state jamming values decreases when the length of the line segment 
/ increases. 
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FIG. 4: The steady state jamming limit against the length of the line segment I for K = 1000. 
The solid line represent the mean field result. 

In Fig. 3 we present the stead-state jamming limit versus the parameter K for I = 2(o) 
and I = 4(D). The solid line is the mean-field prediction 9(oo) = 1 — (1/ ' K) 1 ^ 1 /I for K >> 
1 [18] . Steady-state jamming limits are not coincident with mean-field predictions. Monte 
Carlo results are always smaller than those of mean-field results. Steady-state jamming 
limits increases monotonically when the parameter K increases. 

In Fig. 4 we show steady-state jamming limits against the length I of line segments. 
Jamming limits decrease monotonically when the length I increases. Mean-field results of 
jamming limits (solid lines in Fig. 4) have a minimum point. However, our Monte Carlo 
results do not have a minimum point of the jamming limit. The non-mean field behavior of 
steady-state jamming limits can be understood by the empty site distribution. 

In Fig. 5 we represent the number of empty sites; the total number of empty sites 
N(t), the number of isolated single empty sites N(x0x,t), the number of two consecutive 
empty sites N(x00x, t), the number of isolated empty sites separated by I = 2 occupied 
sites N(x0xx0x, t), and the number of an isolated empty site and two consecutive empty 
sites separated by I = 2 occupied sites N(x0xx00x, t) where x means an occupied site and 
denotes an empty site. The total number of empty sites N(t) = (1 — 0(t))L saturates to 
steady-state limits at long times as shown in Fig. 5 (a) and (c). The early time behavior of 
the empty site is controlled by collective behaviors of the single empty site, double empty 
sites, and higher empty sites. The number of two consecutive empty sites N(x00x, t) shows 
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FIG. 5: The time dependence of empty site distributions with 1 = 2. We present different kinds 
of empty site distributions, the total number of empty sites N (t) (solid line), the number of an 
isolated empty site N(xOx, t) (dotted line), the number of consecutive two isolated empty sites 
N(x00x, t) (dashed line), the number of two isolated empty sites separated by I = 2 occupied sites 
N(x0xx0x, t) (longdashed line), and the number of an isolated empty site and two consecutive 
isolated sites separated by / = 2 occupied sites N(x0xx00x, i)(dotdashed line) (a) for K = 10, (c) 
for K = 1000. Short time behaviors of iV(x0xj;00x, t) (b) for K = 10 and (d) for K = 1000. 

a minimum point around the time t\ and saturates to a steady-state value at t > t\ as shown 
in Fig. 5 (b) and (d). In particular, we observe that N(x0xx0x, t) ^ N(x0x,t)N(x0x,t) 
and N(x0xx00x, t) ^ N(x0x,t)N(x00x,t). The mean-field prediction is based on the ap- 
proximation N(x0xx0x, t) = N(x0x,t) 2 and N(x0xx00x, t) = N(x0x,t)N(x00x,t). We 
conclude that the coverage fraction of the reversible RSA can not explain by the mean-field 
approximation because there is strong collective behaviors of empty site distributions. For 
large f^(for example K = 1000) the total number of empty sites N(t) are controlled by 
N(x0xx0x, t) at the early time t < t\. 

In Fig. 5 (c) we also observe typical four stages of the coverage fraction. At t < t\ the 
coverage fraction increases due to the dominant process of the adsorption. At ti < t < t 2 
the coverage fraction saturates to the jamming limit of the irreversible RSA. In these region 
N(x0x, t) and N(x0xx0x, t) decays slowly up to time t 2 as shown in Fig. 5 (c). At t 2 < t < t% 
the coverage fraction increases as 6{t) ~ log(t) and N(t) ~ (logt) -1 . In this region isolated 
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single empty sites N(xOx, t) and N(xOxxOx, t) decrease logarithmically, while (rrOOx, t) 
decreases very slowly and saturates to a steady-state value. At t > t 3 empty sites distribution 
converges to a limiting value and the empty sites strongly depend on isolated empty sites 
at t > i 3 . 

CONCLUSIONS 

We have observed that the coverage fraction of the reversible RSA does not show the 
mean-field behavior on the one-dimensional lattice. Jamming limits of the coverage frac- 
tion decrease monotonically for the length of the line segment I. The distribution of two 
consecutive empty sites is not production of the distribution of isolated single empty sites. 
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